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of  proposed  theoretical  and  experimental  research,  some  of 
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Cavitation  -  the  formation  and  subsequent  behavior  of 
gaseous  volunes,  usually  in  a  continuous  liquid  phase  -  has  long 
occupied  a  prominent  position  in  engineering  theory  and  practice. 
The  damage  resulting  fron  aquatic  cavitation  say  be  the  preaier 
example  of  the  practical  importance  of  this  phenomenon.  As  a 
result,  the  study  of  naval  architectural  hydrodynaaics  in 
particular  reaains  a  very  active  field  in  which  the  liaited 
reproducibility  and  questionable  direct  applicability  of 
laboratory  work  continues  to  challenge  workers. 

In  reviews  by  Plesset  and  Prosperetti  [65]  and  Acosta  and 
Parkin  [1]  a  natural  concentration  on  aqueous  systens  is 
apparent.  Yet,  even  for  this  specific  Newtonian  fluid,  the  role 
of  geoaetry,  both  nacroscopic  and  aicroscopic,  and  nuaerous  other 
possible  design  variables,  remains  unclear.  Also,  two  competing 
regiaes  which  divide  the  field  exist:  nucleation  or  aicroscopic 
inception,  and  babble  or  cavity  dynaaics.  The  relative 
importance  of  each  regime  varies  according  to  the  specific 
circuastances. 

The  phenomenon  is  further  complicated  by  the  introduction  of 
dilute  aacroaolecular  solutions  into  experiaental  systens. 
Cavitation  inception  in  flow  past  blunt  bodies  is  definitely 
inhibited.  To  aeasure  the  state  of  liquid  flow  at  which 
cavitation  appears,  an  "incipient  cavitation  number,"  e,  has 
been  defined 


dilute 


°i  * 


ps  -  p* 

V2  p  vj 


where  P  is  the  local  free  stream  static  pressure,  P  the  liquid 
s  v 

vapor  pressure,  p  is  the  liquid  density,  and  V  the  free  streaa 

o 

velocity.  Ellis,  £t  al.  [  25,  27]  have  experiaentally  found  a 
reduction  of  by  as  nuch  as  70%  froa  its  value  for  tap  water 
for  a  solution  of  300  ppm  Guar  gum. 

That  this  cavitation  inhibition  is  not  uerely  a  nucleation 
effect  is  demonstrated  by  the  change  in  appearance  of  the 
cavities  upon  fornation.  Ting  [83]  shows  photographs  which 
display  cavities  which  are  saaller  than  their  counterparts  in 
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par*  water.  Hoyt  [43]  and  Ellis,  f_£  al.  [26]  have  also  connented 
on  this  qualitative  change  in  the  natare  of  the  cavities. 

Ellis  and  Ting  [26,  84]  speculated  that  this  cavitation 
suppression  is  a  manifestation  of  the  elastic  properties  of  a 
aon-le vtonian  solution  which  sight  siailarly  inhibit  the  growth 
of  an  individual  cavitation  bubble.  To  test  this  hypothesis, 
they  investigated  the  spherically  sysaetric  growth  of  a  single 
bubble  in  an  otherwise  guiescent  dilute  polyaer  solution  [81,  82, 
84 ].  Agreement  between  theoretical  analyses  and  experiments  was 
good,  but  in  both  the  presence  of  the  nacromolecules  retards 
bubble  growth  only  slightly.  Since  this  retardation  is 
insufficient  to  account  for  the  cavitation  inhibition  in  flows, 
another  approach  aust  be  taken. 

Dilute  polyaer  solutions  can  have  strong,  and  often 
unexpected,  effects  on  the  velocity  field  of  a  flow.  Acoustic 
streaming  is  a  prime  example  [15].  In  this  case  the  direction  of 
the  secondary  flow  generated  by  an  oscillating  rod  is  actually 
reversed  by  the  addition  of  a  snail  anount  of  polyaer. 
Theoretical  analysis  has  shown  that  this  reversal  is  consistent 
with  the  stresses  introduced  by  the  viscoelasticity  of  the 
liquid.  Thus,  there  is  a  strong  possibility  that  the  alteration 
of  overall  flow  field  reduces  cavitation  effects  by  changing  the 
stress  field  around  the  point  of  cavitation,  both  before  and 
after  nucleation.  A  successful  analysis  of  cavitation  inhibition 
is,  therefore,  expected  to  require  inclusion  of  the  overall  flow 
along  with  non-Hewtonian  elastic  effects. 

The  model  system  proposed  here  is  a  single  cavitation  bubble 
growing  in  a  specified  undisturbed  flow  of  a  non-Newtonian  fluid. 
Any  analysis  of  this  systen  involves  a  coupling  between  fluid  and 
flow  which  necessitates  the  synthesis  of  manifold  disciplines. 
One  is  Newtonian  bubble  dynanics,  which  has,  historically,  dealt 
most  completely  with  spherically  syaaetric  flows.  Here,  it  is 
necessary  to  understand  the  relative  importance  of  the  numerous 
physical  effects  and  paraaeters  which  arise.  in  addition,  the 
imposed  flow  will  often  invalidate  any  assumption  of  spherical 
symmetry  so  the  extensive  work  on  small  drop  deformation,  slender 


bod;  theory  and  singularity  analysis  nay  be  applicable.  To 
inclnde  non-Iewtonian  effects  the  various  constitutive 
foraulations  and  formalisations  should  be  evaluated.  Previous 
applications  of  these  constitutive  relations,  particularly  to 
bubbles,  are  especially  important. 

In  subsequent  sections,  each  of  these  subjects  will  be 
treated  in  some  depth.  In  addition,  several  preliminary  results 
attenptinq  to  evaluate  the  iaportance  of  various  physical  effects 
and  to  extend  some  other  results  vill  be  presented.  Further 
steps  of  a  siailar  nature  are  also  proposed.  Finally,  the  sore 
ambitious  plans  for  this  research  are  presented.  These  include 
some  rudimentary  ideas  for  experimental  work  which  will,  it  is 
hoped,  give  guidance  for  theoretical  work  and,  ultiaately, 
confirn  any  analytical  results. 
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II  Uimni  BUBBLE  BIMIICS  and  fi&0£ 


Just  as  most  subsequent  researchers  have  been  motivated  to 
study  bubbles  by  interest  in  cavitation  daaage,  Rayleigh  [72,  50] 
found  this  field  sufficiently  appealing  to  perform  the  first 
analysis  of  bubble  dynanics.  He  solved  the  problea  of  the 
collapse  of  an  empty  cavity  in  a  large  aass  of  liguid. 
■eglecting  surface  tension,  fluid  viscosity,  and  any  thermal 
effects,  and  assuaing  spherical  synaetry,  his  result  can  easily 
be  derived  from  the  equations  of  motion  as 

3  i2  P<R>-  P« 


RR  + 


3  .2 

- R  = 

2 


where  p  is  liquid  density,  p  is  the  pressure  in  the  liquid  far 

m 

from  the  bubble,  and  p(R)  is  the  pressure  at  the  bubble  boundary, 
which  is  also  assumed  to  be  the  pressure  throughout  the  bubble. 
The  specification  of  an  incompressible  fluid  implies  that  the 
velocity  within  the  liquid  is  given  by 

R  *  . 

u(r,t>  =  — rR  (3) 


The  neglected  surface  tension  and  viscosity  terms  can  easily 
included  to  give  a  "generalized  Rayleigh  equation"  [65] 

••  3  .2  1  2a  4|j  . 


RR  ♦ 


3  .2 

- R 
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- [P.  -  P  - 

pi  - 


4jj  . 
— R] 


where  p  is  now  the  pressure  inside  the  bubble  which  is  assumed 
i 

to  be  homogeneous,  and  a  and  p  are  the  surface-tension  constant 
and  liguid  viscosity,  respectively.  In  the  general  situation. 


both  p  and  p  may  be  functions  of  time  (or  of  R)  and  a  and  p  may 

-  i 

not  be  constant.  It  is  also  important  to  eaphasize  that 
spherical  symmetry  is  assumed  to  persist  throughout  the  process. 

For  pure  water,  a  and  p  are  weak  functions  of  other  physical 
parameters,  such  as  tenperature.  This  has  allowed  successful 
modelling  under  the  assumption  of  constant  values  for  those  two 
coefficients.  However,  the  composition  and  behavior  of  p  with 
the  evolution  of  the  bubble  is  not  so  simple.  Plesset  [65]  has 
delineated  two  categories  of  bubble  dynamics  which  he  has 
labelled  gas  bobbles  and  vapor  babbles.  Gas  bubbles  are  those 
cavities  for  which  the  aediua  in  the  interior  is  largely  or 
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completely  a  permanent,  noncondensable  gas.  For  vapor  babbles 
the  gaseous  phase  consists  almost  entirely  of  the  vapor  of  the 
snrronnding  fluid. 

Gas  bubbles  have  been  extensively  investigated,  particularly 
in  the  areas  of  surface  oscillations  and  mass  diffusion  effects 
[see  review  65].  They  have  technological  applications  to  such 
processes  as  plastics  foaming,  but,  for  the  field  of  cavitation, 
these  gas  bubbles  are  much  less  relevant  than  vapor  bubbles. 

The  category  of  vapor  babbles  is  subdivided  into  two  other 
topics.  The  subclassification  is  based  upon  the  extent  to  which 
thermal  effects  alter  p  and  thus,  bubble  behavior.  The  strong 
temperature  dependences  of  equilibrium  vapor  pressure  and  vapor 
density  act  to  significantly  reduce  p  for  a  growing  bubble  as 
evaporation  at  the  bubble  surface  cools  the  interior.  Using  an 
energy  balance,  and  assuming  that  heat  is  supplied  by  a  liquid 
layer  which  has  thickness  comparable  to  the  diffusion  length 
(Dt)  ,  the  "thin  thermal  boundary  layer"  assumption,  this 
temperature  difference  can  be  estimated  as 


„  1 

AT  =  - 

3 

where  L  is  the  latent  heat 
equilibrium  vapor  density  at 


(5) 


H  Py(T)L 
(Dt)*  pC, 

e 

of  evaporation,  p  (T)  is  the 

v 

temperature  T,  D  is  the  thermal 

For 
at 


diffusivity,  and  t  is  the  time  required  to  grow  to  radius  R. 

o  -3  o 

water  at  15  C,  with  B=0.1  cm,  and  t=10  sec,  aT=  0.2  C,  while 
o  o 

100  C,  aT=  13  C.  For  the  former  case,  in  which  the  thermal 
effect  is  not  expected  to  be  important,  the  proper  term  is 
"cavitation”  bubble.  In  the  latter,  and  all  cases  where  thermal 
effects  dominate  inertial  effects,  the  result  is  "boiling"  or 
"vapor"  bubbles. 

The  cavitation  bubble  is  the  simpler  case  since,  with 
complete  neglect  of  thermal  effects,  p  *  p  =  constant.  Then, 
for  a  constant  p  and  neglecting  viscous  effects 

.2  R*  3.2  2  p„-p«.  B.  3  2 o  B0  2 

*  -  <r>  8  4  t  [  i  -  <r>  i  -  -  <r>  ]  (6) 

B  o  3  p  R  pR  B 

where  initial  values  are  signified  by  the  subscript  zero.  For 
Rs  o  and  again  neglecting  surface  tension,  the  Rayleigh  result 
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(7) 


for  tiae  to  complete  collapse  is 

r(5/6)  3w  p  1/2  p  1/2 

t  =  - [ - - - ]  B  *  0.915  ( - - - )  R 

o  r  (1/3)  2(p_-pw)  o  p„-pv  ° 

Quantitative  experimental  confiraation  of  this  value  was  given  by 

Laoterborn  [52].  He  was  able  to  generate  "empty"  cavities  by 

focusing  a  pulse  from  a  Q-switched  ruby  laser  on  a  point  on  the 

interior  of  a  liquid  aass.  Agreement  between  experiment  and 

eguation  (7)  was  excellent,  despite  the  actual  compressibility  of 

the  liquid.  Rote  that  this  compressibility  is  expected  to  become 

significant  in  the  final  stages  of  collapse,  when  the  Hach  number 

of  the  bubble  wall  becomes  large. 

To  analyze  the  more  complicated  vapor  bubbles,  it  is 

necessary  to  solve  the  energy  equation  along  with  the  coupled 

Rayleigh  equation.  Plesset  and  Zwick  [66,  67]'  found  that  under 

the  assumption  of  a  thin  thermal  boundary  layer,  the  temperature 

at  the  bubble  wall  is  approximately 

D  1/2  Lp'  t  8l(x)  R  (x) 

T  »  T  -  ( - )  — dx  (8) 

-  *  k  x  =  0  [/*  R4<y)  dy]*fc 

where  T  is  the  temperature  at  infinity,  D  is  the  liquid  thermal 

am 

diffusivity,  L  is  the  latent  heat  of  evaporation,  p*  is  vapor 
density,  and  k  is  the  thermal  conductivity  of  the  liquid.  The 
growth  of  the  vapor  bubble  is  governed  by  numerous  parameters, 
including  A  which  gives  a  linear  approximation  to  the 
relationship  between  vapor  pressure  p  (T)  and  T 

Pv(T)  -  P.  V 

- —  =  A  (T  -  T  )  (9) 

p  b 


where 


is  the  boiling  temperature  at  pressure  P  .  Other 


constants  are  defined  as 


a  -  {' 


(70  A) 


3  kR0a 


D  1/2 
( - ) 

■a 


(10B) 


mnd  $,  which  is  implicitly  given  as  a  root  of 

2  1/2 

$  ♦  3p  («£)  -1=0  (10C) 

If  the  mechanism  by  which  growth  is  triggered  from  an  initial 
equilibrium  radius  B  is  a  constant  heat  source  per  unit  volume 


of  strength  a,  such  that 


u  =  at  (1 OD) 

then  the  Rubble  evolution  can  divided  into  four  regiaes  [93]. 
The  two  most  easily  understood  are  the  first,  a  "delay"  period, 
and  the  fourth,  where  growth  approaches  an  asyaptotic  Unit. 

The  initial  growth  occurs  when  radial  increase  is  so  rapid 
that  the  change  in  bulk  temperature,  which  follows  from  equation 
(10D) ,  is  insignificant.  To  first  order,  the  solutions  to  the 
governing  eguations  are 

.  2r 

R  =  H  (1  ♦  exp  {apt)  ]  (11A) 

R 5a*  2rO-0  ) 

T  (F)  =  T - . . . -  exp  (aftt)  (11B) 

-  A  ft  ^3ft%*■^) 


where 


3  2 

y  =  (ADa)  / (a  R  k)  (12) 

o 

The  delay  period  is  very  short,  but  results  in  a  new  value  of  R, 

slightly  larger  than  R  ,  froa  which  "iaportant"  growth  begins. 

o 

Bhen  R>>R  ,  the  heat  diffusion  has  a  limiting  effect  on  R. 
o 

A  consistent  solution  to  the  equations  of  motion  and  energy  which 
takes  into  account  the  delay  period  is  generated  using  an 
asymptotic  expansion.  The  leading  terms  of  this  solution  are 


2  «t  1/2  -1/2 

B  =  B  ( - )  <-— )  [  1  ♦  0(t  )  ] 

O  wp  3 


T  =  T  - 
o 


a  K  _ 


-1/2 

[  1  ♦  0(t  )  ] 


(13A) 


(13E) 


The  two  intervening  regiaes  can  be  thought  of  as  interaediate 
behaviors  necessary  for  patching  together  the  first  and  last  tiae 
periods.  Hotice  that  Rat  is  a  decreasing  function  of  tiae 
in  the  asymptotic  regime. 
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FIG. 1  [231 


The  large  aagnitude  of 

the  thermal  effect  on  boiling 

babbles  is  illustrated  in 

figure  1.  Here  the 

theoretically  predicted  radius 

vs.  time  profiles  for  the 

Plesset  analysis  is  shown 

along  with  the  Rayleigh  result 
o 

for  water  at  103  C.  Excellent 
experimental  agreement  with 
the  Plesset  prediction  was 
found  by  Dergarebedian  [23]. 
He  used  a  Q-switched  ruby 
laser  to  generate  bubbles 
under  conditions  closely 
approxinating  the  theoretical 
situation. 


Zwi ck  and  Plesset  [93]  have  also  analyzed  the  collapse  of  a 
bubble  including  the  heating  effect  caused  by  condensation  of  the 
vapor. 
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Here  numerical  techniques 
are  definitely  required 
since  the  expected  large 
temperature  differences 
do  not  allow  simple 
analytic  expressions  for 
vapor  pressure  or  vapor 
density,  evidently,  the 
time  scale  for  collapse 
is  so  snail  that  the 
dynamics  of  the  process 
vary  only  slightly  fron 
the  Rayleigh  solution 
(See  figure  2)  . 


FIGURE  2  [93,  P.3238] 

These  analyses  all  contain  the  premise  of  a  spherically 
symmetric  bubble.  Plesset  [61  ]  initially  performed  a  linear 
stability  analysis  on  the  nearly  spherical  interface  between  two 
immiscible,  incompressible,  inviscid  fluids.  He  assumed  an 
initial  drop  shape  given  by 

r  =  R  ♦  a(t)  Y  (ih) 
s  n 

where  Y  is  a  spherical  harmonic  of  degree  n  [50,  5tt,  88,  92]  and 
a  is  initially  small.  The  analysis  proceeded  conventionally  in 
order  to  determine  the  conditions  for  which  a(t)  will  grow, 
implying  the  shape  is  unstable,  and  those  for  which  it  will  not, 
implying  stability,  one  finds  that 


a(t)  *  «<t)  (Ji£-)3/2 
»ith  a  -  G(t)  «*0 

•  t  M 

3  R  R  3 

t  r*-  *  + 


(15) 

(16A) 


G(t) 


n  (n-1)Pt.-<n  +  1)  (n*2)p, 
npj  ♦  (n*1) 


] 


(16B) 


(n- 1 )  n  (n+  1 )  (n*  2)  */R  * 

R[n ♦  (n+1)^] 

where  a  is  the  surface  tension  constant  and  p  and  p  are  the 

1  2 

fluid  densities  of  the  interior  and  exterior  liquids, 
respectively.  Thus 

6 (t)  <0  promotes  stability 

G  (t)  >0  proaotes  instability. 

This  aeans  that  surface  tension  always  has  a  stabilizing  effect, 
and  from  (15)  ,  increasing  R  or  bubble  growth  is  seen  to  have  a 
stabilizing  effect,  while  the  opposite,  bubble  collapse,  promotes 
instability. 

Plesset  and  Hitchell  [64]  then  performed  a  more  involved 
analysis  for  a  vapor  cavity,  neglecting  the  density  of  the  vapor 


and  the  viscosity  of  the  vapor  and  the  liquid.  Results  showed 

that,  for  an  expanding  vapor  cavity,  if  |a(0) |/R  <<  1,  then 

o 

|a(t)|/R(t)  «  1.  For  a  collapsing  cavity  distortion  amplitudes 

remained  small  so  long  as  1.0  >  R/R  >  0.2,  but  as  R-»0,  a(t) 
-1/4  o 

increased  as  R  .  Thus  the  spherical  shape  is  unstable  for  the 
later  stages  of  bubble  collapse.  These  trends  are  valid  even 
when  a  small  viscous  effect  is  included  in  the  treataent  [70]. 
Viscosity  does  tend  to  damp  the  growth  in  amplitude  of  the  higher 
order  harmonics. 

To  determine  bubble  shape  stability,  no  characteristics  of 
the  deformation  beyond  its  existence  and  magnitude  were 
necessary,  a  direct  result  of  the  far-field  conditions.  The 
bubbles  were  assumed  to  exist  in  an  infinite  fluid,  quiescent 
apart  from  the  direct  effects  of  the  cavity;  this  is  a 
spherically  symmetric  geometry  giving  an  equilibrium  shape  with 
similar  character.  Of  course,  this  is  not  the  only  circumstance 
in  which  a  bubble  may  arise.  There  may  be  conditions  imposed  on 
the  fluid  away  from  the  bubble,  such  as  an  elongational  flow  or  a 
solid  boundary,  or  the  bubble  may  be  initialized  as  nonspherical 
and  a  detailed  description  of  its  evolution  desired.  Bubbles 
have  provoked  some  study  in  this  regard,  but  aore  often  liquid 
drops  and  solid  bodies  have  been  exaained,  giving  rise  to  general 
techniques  for  nonspherical  shapes. 


For  babbles  end  all  initially  spherical  bodies  in  an  iaposed 
flow,  the  primary  method  of  study  has  been  expansion  of  the 
velocity  profile  in  terms  of  spherical  harmonics  [76].  Ose  of 
this  aethod  is  predicated  upon  an  assuaption  of  creeping  flow,  at 
least  in  the  immediate  vicinity  of  the  babble.  Onder  this 
restriction,  Happel  and  Brenner  [38,  sect  3-2]  offer  a  practical 
presentation  of  the  use  of  spherical  harnonics,  Laab's  General 
Solution  [50].  The  pressure  field  nast  satisfy 
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0 


This  suggests  expanding  p  as  a  series  of 

harnonics  p  . 

n  • 


p  =  £ 

n*-« 


n 


Then  Laab's  General  Solution  for  velocity  y  is 

■  (n+3)  2 

v  -  r  [7x(rX  )  ♦  7e  ♦  - - - — -  r  7p 

n=-  •  ~  n  n  2*i(n*  1)  ( 2n+  3)  n 


(18) 

solid  spherical 
(19) 


n 

-  rp  ] 

v(n+1)  (2n+3)  n 

(20) 


where  x  and  *  are  each  solid  spherical  harnonics  given  as  the 
n  n 

series  of  solutions  to 

2 

7  y*  *  0  (21 A) 

V.y  =  0  (2 1 B) 

Happel  and  Brenner  proceed  to  describe  the  determination  of  the 
p  ,  X  and  ♦  for  given  boundary  conditions. 

Einstein  [24]  was  able  to  calculate  an  equivalent  viscosity 
Tj*  for  a  dilute  suspension  of  rigid  spheres  without  using 
spherical  harnonics.  His  result 


5  3 

tj*  =  »j  (1  ♦  —a  )  (22A) 

2 
3 

a  <<  1  (22B) 

where  is  the  viscosity  of  the  dispersing  fluid  and  a  is  the 
volume  fraction  of  spheres,  was  rederived  using  a  harmonic 
expansion  by  Frohlich  and  Sack  [32].  They  were  also  able  to 
calculate  the  viscous  and  elastic  properties  of  a  substance 
equivalent  to  a  dilute  suspension  of  elastic  spheres  in  a  steady 
elongational  flow.  The  non- rigidity  of  these  elastic  spheres 
introduces  the  need  to  calculate  the  defornation  which  would 
result  if  these  same  stresses  were  present  for  the  elastic 
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spheres.  In  essoining  this  to  give  an  accurate  description  of  the 
non-sphericity,  the  defornation  is  required  to  be  small  and  not 
to  pertorb  the  flow  significantly.  In  considering  the  behavior 
of  a  viscous  drop  in  the  flow  of  a  viscous,  lewtonian  fluid,  two 
important  parameters  are  found  to  be 

X  *  m*/*»  (2 3A) 

o 

k  *  ofti  Ga  (23 B) 

o 

where  **  is  the  the  viscosity  inside  the  drop,  u  *  that  of  the 

surrounding  fluid,  a  the  surface-tension  constant,  a  is  the 

radius  of  the  drop  and  G  some  measure  of  flow  strength.  Taylor 

[76]  discussed  the  result  for  a  siaple  shear  flow,  and  found  two 

limiting  cases.  For  the  first,  with  A=0(1)  and  k>>1,  terned  the 

interfacial  tension  doninated  case,  the  drop  deforas,  to 

o 

order  (1/k) ,  into  a  spheroid  with  its  najor  axis  45  to  the  flow. 
For  the  second,  with  k=0(1)  and  A»1,  to  0(1/*)#  the  major  axis 
is  aligned  with  the  flow.  Taylor  and  icrivos  [79]  found  that  a 
falling  drop  was  deforaed  into  an  oblate  spheroid  when  surface 
tension  dominates,  and  has  spherical  caps  as  surface  tension 
becomes  less  important.  They  also  performed  higher  order 
calculations  by  re-evaluating  the  flow  around  the  deforaed  drop. 

A  general  tiae  dependent  flow  containing  a  viscous  drop,  and 
other,  more  specific,  flow  situations  have  been  analyzed  by  Cox 
[20],  Frankel  and  Acrivos  [31],  Barthes-Biesel  and  Acrivos  [4], 
and  others  (  36,  37].  The  viscosity  ratio  k  (23a)  and  ratio  or 

surface  tension  to  flow  strength  k  (23b)  reaain  iaportant 
parameters.  In  the  notation  of  Cox  [20],  the  deformed  surface  is 
given  by 

r  *  1  ♦  vf(r  /r)  (24) 

i 


where  r  is  a  small  parameter.  The  function  f  is  expanded  in 
teras  of  spherical  surface  haraonics  as 

9l 


f  *  F 

pq 


9rfdt% 


(1/r)  ] 


r*1 


(25) 


Bakimi  [36]  presents  experimental  results  which  compare  well  with 
theoretically  derived  values  for  F  .  His  work  aade  use  of  an 

pq 


2 


orthogonal  rhooaotor  to  gonorata  flow  fields  with  variable 
contributions  from  vorticity  and  deformation.  Despite  this 
success,  there  is  "sose  confusion"  anong  various  workers  about 
the  agnations  to  be  solved  in  deteraining  the  P  .  Ballison  [71] 
atteapts  to  unify  these  differing  results,  lid  in  so  doing 
concludes  that,  nithin  the  range  of  validity  of  this  analysis, 
there  are  at  least  three  special  cases  in  which  an  equilibrium 
shape  nay  be  achieved.  These  special  cases  are: 
leak  flow  -  k»1  with  A«o(k) 

High  viscosity  drops  -  a>> 1  with  k*o(x) 

High  viscosity  with  comparably  weak  flow  -  k,  a>>  1 
For  bubble  dynamics,  only  the  first  case  may  be  useful. 


o 


Comparison  of  experimentally  determined  bubble  dupM  (open  circles)  on. 
ootlapeo  of  a  epherioal  bubble  near  a  plana  aoli  J  wall  with  theoretical  curve*  taken  from 
Pleeeet  a  Chapman  (1971)  (eolid  curves).  The  framing  rata  ia  >00 000  fiamee/e,  the  maxi- 

num  bubble  radius  R _ m  i  t  mm,  the  distance  of  the  bubble  centra  from  the  wall 

b  e  S-B  mm  and  bIR _ e  1-6. 

Curve  01  S  346*7 

Time.  0  0-726  0-126  0-961  0-991  i-010  1-02$  1-03C 

(p,  is  the  density  of  the  liquid  and  Ap  ia  the  aonatant  difference  between  the  ambient 
liquid  preasure  and  the  pressure  in  the  cavity.) 

FIGURE  3  [63,  P.396] 
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The  presence  of  a  solid  boundary  in  the  vicinity  of  a  babble 
alters  the  flow  created  by  the  cavity  dynamics.  In  the  collapse 
process  a  cumulative  jet  is  formed  which  may  be  one  of  the  major 
causes  of  the  destructive  action  of  cavitation.  Bsperiaentally, 
the  jet  has  been  observed  in  laser-induced  bubbles  [53],  with 
qualitative  agreement  achieved  with  numerical  calculations  [65] 
(See  figure  3) .  The  theoretical  aodel,  however,  does  not 
generate  jet  speeds  sufficient  for  the  destruction  encountered  in 
practice,  so  that  other  factors  such  as  nonspherical  initial 
shape  have  been  investigated  [70]. 

An  initial  nonsphericity  is  important  even  in  the  absence  of 
an  iaposed  flow  or  solid  boundary. 

The  requisite 
aatheaatics  has  been 
put  forth  by  Hsieh 
[45],  who  includes 
equations  for  heat 
and  aass  transfer 
effects, 

coapressibility  and 
suggests  variational 
aethods.  Simplified 
numerical  results  by 
other  authors  [17] 

Bubble  ivrfacee  fram  cate  A 

show  that  an 
initially  prolate 
spheroid  (figure  4a) 
will  form  two  jets 
upon  collapse,  while 
an  oblate  spheroid 
(figure  4b)  leads  to 
a  dumbbell  form. 

These  nonlinear 

results  were  I  ■•ww# •*»*■*••  •** a 

compared  with  the  1  FIGURE  4  ri71 

linear  theory 
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developed 


by  Plesset  eed  Bitcbell  [64]  la  their  detailed  stability 
aaalysis  and  good  agreeaent  eas  found  until  the  final  stages  of 
collapse. 

The  success  of  the  linear  analysis  in  predicting  bubble 
behavior  even  for  large  defornations  is  a  pleasant  surprise. 
Bowever,  such  analyses,  which  begin  from  perfect  spheres,  cannot 
be  expected  to  have  such  a  large  range  of  validity.  Even  the 
cited  linear  analysis  eventually  becoaes  invalid.  Nonetheless, 
these  shapes  are  important;  there  is  a  strong  possibility  that 
bubble  dynamics  will  generate  highly  deformed  shapes  in  certain 
isposed  flows  since  the  snail  vicosity  ratios  a  (23a)  present  do 
not  result  in  a  "si all  detonation"  eguilibriun,  in  the  sense  of 
Ballison  [71],  unless  the  flow  strength  is  very  low,  or  surface 
tension  dominates.  The  most  common  method  for  dealing  with 
highly  deforaed  bodies  is  called  slender  body  analysis.  The 
theory  neglects  inertial  effects  and,  as  the  naae  suggests,  has 
been  created  for  flows  around  long,  slender,  solid  bodies 
characterized  by  snail  thickness  ratios  «. 


«c  *  b/1  «  1  (26) 
where  the  cross-sectional  radius  is  of  order  b,  and  1  is  the 
length  of  the  body. 

Slender  body  theory  was  initiated  by  Burgers  [13],  when  he 
attempted  to  determine  the  force. on  a  long,  slender  ellipsoid  ot 
revolution  at  rest  in  a  uniforo>  translational  motion  0.  He  began 
the  analysis  by  assuaing  the  disturbance  produced  by  the  body  was 
like  that  which  would  result  fron  a  line  of  force  on  the  synnetry 
axis,  which  is  of  magnitude 


2  4 

f  (z)  =  A  ♦  A  (z/a)  ♦  A  (z/a)  |z|<a 

o  2  4 

f(z)  -  0  |z| >a 


(27) 


where  a  is  the  semi-major  axis,  z  is  the  distance  along 

from  the  center  of  the  ellipsoid,  and  A  ,  A  ,  and 

0  2 

constants.  He  was  able  to  derive 

a  4«ua0 _ 

ln(2a/b)  -  0.5 


the  axis 

A  are 
4 


(28) 
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which  is  the  Uniting  form  of  the  re  salt  obtained  for  spheroids 
of  finite  cross-section.  Subsequent  efforts  in  this  field  by 
other  workers  have  extended  this  type  of  analysis  to  non "Uniform 
cross-sections  [6],  shear  flows  [21(2)],  and  even  curved  axes 
[21(1)  ].  In  all  cases,  the  intended  result  is  the  function  F(z) 
which  best  matches  boundary  conditions,  and  is  often  obtained 
using  matched  inner  and  outer  expansions.  From  F(z),  the  forces 
and  torques  on  the  body  can  be  calculated. 

k  sinqle  function  F(z)  can  give  exact  results  only  in  the 
linit  as  *~>0.  This  function  F(z)  can  equivalently  be  considered 
the  distribution  along  the  axis  of  point  singularities  in  the 


flow  field,  or  Stokeslets  [5].  The  velocity  and  pressure  field 
due  to  a  Stokeslet  P  [11]  is 

j  . 

U.  =  [— *>  ♦  — /  ]  (29 k) 

i  Bmu  r  r* 


(29B) 


For  bodies  of  finite  cross-section,  it  is  only  possible  to 
generate  exact  solutions  to  Stokes  flow  problems  by  determining 
the  proper  spatial  distribution  of  Stokeslets  &nd  the  other 
fundaaental  singularities  (figure  5) ,  which  can  be  derived  fron 
the  derivatives  of  the  Stokeslet  [11].  Thus,  just  as  snail 
deformation  methods  may  be  valid  for  finite  deformations,  the 


slender  body  singularity  technique  is  extendable  to  finite 
thickness  ratios.  Host  shapes  are  analyzable,  at  least  for 
Stokes  flow  and  rigid  bodies. 


XTOKSS  OOUWLZT  STRCSSLCT  WOTUET 

_ ( COUHt  T) 

FIGURE  5  [  11,  P.25] 


16 


Ill  BOW— BBHTOBI AH  EFFECTS 

Many  instances  of  seemingly  anomalous  behavior  of 
non-Hewtonian  fluids  are  well-docunented.  Two  of  the  aost 
commonly  cited  examples  are  rod  climbing  or  the  Reissenberg 
effect  [8,  P.92],  and  extrudate  or  "die**  swell  [8,  p.102]. 

Schowalter  [74]  emphasizes  the  "counter-intuitive**  nature  of 
these  effects,  and  explains  that  this  is  a  conseguence  of  the 
Bewtonian  character  of  that  very  cosson  fluid,  water,  which  is 
the  basis  for  such  intuition.  while  water,  and  all  truly. 
Bewtonian  fluids,  exhibit  no  elastic  forces,  non-Newtonian  fluids 
often  do,  and  the  existence  of  such  forces  can  be  used  to  explain 
aany  of  the  unexpected  effects. 

Several  of  these  effects  can  be  viewed  as  exhibitions  of 
fluid  properties  which  are  also  relevant  to  cavitation 
inhibition,  particularly  those  properties  which  involve  a  free 
surface  or  result  in  a  macroscopic  alteration  of  the  flow  and/or 
pressure  field.  In  the  first  category,  exaaples  include: 

jet  stability  -  in  a  Newtonian  fluid,  breakup  of  a  free 

jet  occurs  in  clearly  defined  waves.  Such  definition 
is  not  apparent  for  fluids  which  exhibit  elastic 
properties  [74,  P.3]. 

open  channel  flow  -  for  many  non-Newtonian  fluids,  the  free 
surface  of  the  liguid  when  flowing  down  an  open, 
inclined  channel  is  convex,  it  is  virtually  flat  for 
the  Bewtonian  case  [8,  P.105]. 

Demonstrations  of  the  second  type  of  effect  are  necessarily  sore 
involved,  reguiring  some  flow  visualization  or  measurement 
technigues.  Tet  there  are  good  exaaples,  including: 

secondary  flow  induced  by  an  oscillating  cylinder  -  as 

previously  cited,  Chang  [15]  observed  reversal  of 
secondary  flow  (figure  6)  for  polyacrylaaide  solutions 
as  dilute  as  20-30  ppn.  Re  also  found  the  situation 
amenable  to  analysis. 
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Secondary  flow  produced  by  a  long  cylinder  oscillating  normal  to  its  axis.  The 
cylinder  is  viewed  on  end  and  the  direction  of  oscillation  is  shown  by  the  double  arrow.  (§)  A  water/ 
glycerin  mixture  moves  away  from  the  cylinder  along  the  axis  of  oscillation.  (P)  The  direction  of  the  secon¬ 
dary  flow  is  reversed  when  100  ppm  polyacrylamide  (Separan  AP  30)  is  added  to  the  water/glycerin 
mixture. 


FI<Vld*  *S  [16] 


Particle  petition  vs.  time  for  Newtonian  fluid  (fluid  4).  Portide  positions  vs.  time  for  viscoelostic  fluid  (fluid  3). 


flow  past  a  flat  plate  - 
experimental  results 
{41]  show  definite 
changes  in  particle 
position  (figure  7)  and 
fluid  velocity  (figure 
8)  when  a  Newtonian 
liguid  is  replaced  by  a 
viscoelastic  fluid  in  a 
flow  with  this  siaple 
geometry. 


- viscoelastic 

- NEWTONIAN 

»4*W»I 


0 

DISTANCE  OOWN  THE  PLATE  <X) 


Schematic  diagram  of  o burred  kiae- 
malic  kekariar  for  eiicoelactic  fluids 


FIGDEE  8  [41,  P.257] 


Finally,  a  aost  puzzling  instance  of  non-Newtonian  behavior,  and 
one  which  actually  involves  bubbles,  arises  in  the  recent 
experimental  work  of  Hassager  [40].  He  photographically 
deaonstrates  the  existence  of  a  knife-edged,  asyaaetric  "negative 
wake"  behind  a  bubble  rising  in  a  non-Newtonian  fluid.  There  are 
aany  unanswered  questions  here,  aost  of  which  have  relevance  to 
cavitation  inhibition. 

To  theoretically  analyze  the  flow  of  any  aaterial,  it  is 
necessary  to  construct  a  model,  which,  if  successful,  will 
"accurately"  predict  the  stress  field  in  that  material  or  class 
of  materials  under  varyinq  conditions  of  notion.  Such  a  aodel  is 
called  a  constitutive  equation  and  is  often  formulated  to 
determine  the  stress  tensor  £  for  given  material  conditions, 
present  and  past.  The  constitutive  equation  must  satisfy 
numerous  physical  principles  and  ideally  will  cover  a  broad  range 
of  materials,  but  this  breadth  should  not  iapair  its  practical 
applicability. 

For  non-Newtonian  fluids,  any  attempt  to  formulate  a 
generally  applicable  constitutive  equation  is  coaplicated,  and 
often  stymied,  by  the  negative  nature  of  the  definition  of  this 
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class  of  satter  - 

■oa-lewtoniap  Fluid:  ill  fluids  which  are  not  Newtonian. 
Under  the  constraint  of  incoapressibility , 

?.v  =  0  (30) 

where  v  is  the  vector  velocity  in  the  fluid.  The  constitutive 
equation  for  a  Newtonian  fluid  is 

S  *  -pi  *  T  (31  A) 

*  “Pi  <tv)T]  (3 1 B) 

where  p  *  -(1/3)  tr  S,  T  is  the  extra  stress  tensor,  and  ji  is  the 
Newtonian  fluid  viscosity.  For  a  sufficiently  sinple  flow,  the 
essential  feature  of  this  relation  can  be  rewritten  as 

r  *  ne  (3 2 A ) 

where  r  is  now  the  extra  stress,  e  is  strain,  and  the  dot  denotes 
time  differentiation.  A  similar  relationship  for  a  purely 
elastic  aaterial  is 

r  =  Ke  ( 32B) 

where  K  is  the  elastic  nodulus.  The  presence  of  viscous  forces 
and  elastic  recovery  in  aany  non-Newtonian  effects  suggests  the 
combination  of  these  two  constitutive  relations  to  describe 
non-Newtonian  behavior,  connecting  elements  exhibiting  viscous 

damping  and  elastic  recovery  in  series,  the  Haxwell  element, 

figure  9,  arises,  with  its  constitutive  relation 

r  ♦  At  =  u  e  (33) 

o 

Other  simple  combinations  of  the  elastic  and  viscous  elements 
result  in  other  linear  viscoelastic  equations,  including 

r  ♦  At  =  u  (e  ♦  A'e)  (34) 

Further  sophistication  is  possible  by  including  multiple  values 
of  viscosity  and  aodulus;  the  adjustaent  of  these  values  allows 
better  agreeaent  with  certain  experimental  data. 


« — QF — wv — • 


FIGURE  9  Haxwell  Bleaent 
The  two  equations  (33)  and  (34) 
representations.  The  Haxwell  eleaent  can  also 
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are  differential 
be  described  in 


integral  fora  as 


r(t)  «  2  (m/a)  £  exp[-  (t-t  •)  /a]  e(t*)  dt*  (35) 

This  integral  form  introduces  the  notion  of  a  memory  faaction. 
The  exponential  kernel  of  the  integrand  weights  the  contributions 
to  r  (t)  by  the  strain  e(t')  at  past  tiaes  t*.  By  its 
mathematical  properties,  this  exponential  memory  function  places 
aore  eaphasis  on  recent  values  of  a.  This  distribution  satisfies 
a  common-sense  concept  formally  stated  as  the  principle  of  fading 
neaory . 

These  linear  viscoelastic  models,  which  eaploy  siaple  tiae 
differentiation,  cannot  be  generally  valid  because  they  do  not 
satisfy  another  prescription  for  constitutive  equations,  the 
principle  of  material  objectivity.  The  mathematical  criteria  for 
satisfaction  of  this  principle  are  complex,  but  the  concept  is 
siaple;  the  behavior  of  a  aaterial  as  calculated  for  a  given 
aodel  aust  be  independent  of  observer  or  fraae  of  reference.  The 
lack  of  "frame- indifference"  of  eguation  (33)  means  that  the 
calculated  response  of  a  material  which  obeys  this  equation,  when 
placed  on  a  rotating  table,  would  vary  with  the  rate  of  this 
rotation,  even  when  inertial  and  relativistic  effects  are 
neglected. 

Another  way  to  view  this  difficulty  with  an  equation  like 
that  for  the  Maxwell  model  is  to  recognize  the  ambiguous  nature 
of  any  measure  of  relative  strain  for  non-infinitesimal 
deformation.  To  illustrate,  consider  a  bar  of  length  1  which  is 
stretched  to  length  1*a1,  figure  10. 


FIGURE  10 
STRAIN 


In  the  Unit  of  inf initesiwal  strain,  as  al-hO,  the  relative 
strain  is  given  by  two  equivalent  expressions 

Al  Al 


However,  if  a1  increases,  these  two  expressions  are  no  longer 
equivalent,  and  there  is  no  natural  choice  between  them  which  is 
generally  valid.  To  use  such  measures  without  discrimination 
implicitly  restricts  their  validity  in  flow  to  inf initesiaal 
relative  strain. 

The  restriction  to  inf initesiaal  relative  strain  is  a  first 
example  of  the  specificity  of  certain  models  to  limited  types  of 
flow.  The  intiaacy  of  the  relationship  between  flow  and  material 
behavior  is  illustrated  very  well  by  the  simplification  of  § 
which  obtains  when  the  tensor  is  evaluated  for  a  coaaon  class  of 
flows,  steady  laminar  shear  flows.  Steady  laminar  shear  flows 
are  characterized  by  a  geoaetry  which  grants  a  natural  orthogonal 
coordinate  system.  The  flow  is  restricted  to  having  a  velocity 
in  only  one  coordinate  direction,  labelled  1,  which  varies  with 
respect  to  only  one  other  direction,  2.  If  the  extra  stress  J  is 
defined  for  an  incompressible  non-Newtonian  fluid  as  it  was  for 
the  Newtonian  case,  (31A) 

?(t)  =  S(t)  ♦  pi  ( 37  A) 
1 

where  p  =  -  tr  S  (37B) 
Then  restrictions  on  T  can  be  found  from  components  of  the 
general  equation  of  notion 


Dv 

♦  div  T 

p  —  = 

p9.  ”  grad  p 

(30) 

Dt 

The  result  is 

1  till)  t(12) 

0  1 

T  = 

t(21  t  (22) 

0 

|  0  0 

t  (33) 

(39  A) 

tr  T  = 

0 

(39B) 

Also,  T  can  only 

vary  with 

the  shear  rate  «  and  is 

symmetric. 

Thus,  three  material  functions,  r,  N  ,  the  first  normal  stress 

difference  and  ■  ,  the  second  noraal  stress  difference, 

2 

completely  specify  the  tensor  T.  In  one  notation  [74],  they  are 
given  by 

t  («)  =  t  ( 12)  =  t  (21)  (40A) 

N^U)  =  t(11)  -  t  (22)  (40B) 

N2<«)  =  t  (22 )  -  t  (33) 


(40C) 


k  ®  dv  /dx 
1  2 

For  a  Newtonian  fluid,  both  noraal  stress  differences 


(UOD) 
are  zero. 


It  is  the  non-zero  value  of  these  quantities  which  is  manifest  in 
such  instances  of  "counter-intuitive"  behavior  as  the  Veissenberg 
effect.  These  values  can  be  seasured.  coaaonly  eaployed 
viscometer  flows  are  analyzable  as  steady  laminar  shear  flows, 
and  have  predictive  value. 

The  characterization  of  a  non-Newtonian  fluid  afforded  by 
aeasureaent  of  the  three  aaterial  functions  (40),  although  useful 
and  relatively  siaple,  is  not  coaplete;  the  degeneracy  in  T  and 
exclusive  dependence  on  shear  rate  *  is  flow-specific.  The 
simple  linear  viscoelastic  aodels  are  also  useful,  but 
unsatisfactory  in  a  general  situation.  These  complications 
suggest  a  "retreat"  to  more  basic  properties  and  to  a  foundation 


for  a  constitutive  equation  or  nodel.  Noll  [55,  85,  86]  has 

presented  the  aost  elegant  treataent.  His  foraulation  for  the 

most  general  constitutive  equation  for  a  body  is 

S(X,t)  =  5  (X)  (41) 

“  t,X 

where  Sr  (X)  is  a  functional  of  the  aotion  of  the  body  iron  all 
t,X 

past  tiae  to  the  present  tine  t.  By  invoking  physically 

reasonable  assumptions,  and  retaining  the  defining  equations 

(37a)  and  (37b) ,  Noll  postulates  the  constitutive  relation  for  an 

incoapressible  siaple  fluid 

T(t)  =£■  [C  (s)  ]  (42) 

s=0  (t) 

where  is  a  functional,  £  (s)  is  related  to  the  relative 

-  (t) 

deformation  qradient,  and  s  is  a  dummy  variable  used  to  denote 

the  dependence  on  all  past  tiae.  This  constitutive  eguation  is 

simple  in  form,  elegant,  and  contains  only  a  few  restrictions. 

If  more  assumptions  are  made,  eguation  (42)  assumes  a  more 

specific  fora.  Bivlin  and  Ericksen  [73]  restricted  the  tiae 

dependence  of  T  to  s=0  (near  the  present)  and  performed  a  series 

expansion.  In  their  notation 

I(t)  -  1(1  .  A,...,*  )  (43) 

1  2  n 

where  the  1  are  the  Bivlin-Ericksen  tensors,  functions  of  the 
i 
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(44) 


velocity  gradient  £,  given  by 


n 


T  n- 1  n  T 

It  ♦  St  *  .  E  (.)  L.  L 
n  n  i= 1  i  i  (n-i) 


A  fore  such  as  (43)  ,  though  sore  restricted  than  that  of  (42)  is 
still  not  very  useful,  and  special  cases  aust  be  examined.  If 
the  dependence  on  all  Rivlin-Ericksen  tensors  higher  than  the 
second  is  suppressed,  and  F  Beets  certain  eatheeatical 
conditions,  it  can  be  expanded  as 


2  2 

T  =  a  ♦  a  i  ♦  «  4  *  a  A  ♦  at  (A  &  ♦  h  h  )  * 

0  1  1  2  2  1  3  2  4  1  2  2  1 


(45) 


Despite  all  of  these  restrictions,  equation  (45)  is  still 
cumbersome  and  in  practice  it  is  often  truncated  after  two  terns. 
This  yields  the  equation  for  a  "second  order  fluid,"  which  has 
been  found  to  be  of  value  in  sone  circumstances. 

Another  approach  to  constitutive  equations  has  been 
notivated  by  the  usefulness  of  the  linear  viscoelastic  equations 
and  the  existence  of  soae  physical  weaning  for  the  constants. 
For  general  situations,  it  is  necessary  to  modify  these  equations 
so  that  they  satisfy  material  objectivity.  Oldroyd  [56] 

recognized  that  the  problee  originated  in  the  taking  of  the  time 
derivative  for  a  coordinate  basis  which  is  fixed  in  the 
"laboratory"  reference  frame.  He  suggested  manipulations  with 
respect  to  coordinates  inherent  to  the  body.  The  result  has  been 
the  formulation  and  use  of  many  specialized  types  of  derivatives 
which  are  ieplicitly  naterially  objective.  These  include  the 
codeforeational  derivative  and  the  corotational  or  Jaumann 
derivative.  The  Jaumann  derivative  is  given,  in  tensor  index 
notation,  by 

&  3 

XVij 

where  v  is  the  velocity  and  w  vorticity.  Hany  constitutive 
equations  have  been  formulated  using  these  derivatives  [8], 
including  the  corotational  Jeffreys  model 

i  ♦  -  *-  «  ♦  vb-b 

where  g  is  the  rate  of  deforest  ion.  From  consideration  of  the 


k  k  k 

—  a  ♦  v  a  -  w  a  -  u  a 

3t  ij  ij,k  i  kj  j  ik 


(46) 


(47) 
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applicable  mathematical  properties,  Oldroyd  added  nonlinear  terms 
in  T  and  £  to  give  the  Oldroyd  eight-constant  aodel  which,  in 
Cartesian  coordinates,  is 


(1  ♦  h 


— )  t  ♦lit  d  —  n(t  d  ♦  t  d  )  ♦  v  t  d  8 
t  ij  o  kk  ij  1  ik  kj  kj  ik  1  kl  kl  ij 


& 

*  2  |i  [  d  ♦  a  v  d  -  2*,  d  d  ♦ 
1  ij  2  &t  ij  2  ik  kj 


ii  d  d  s  1 
2  kl  kl  i j 1 


(4  8) 


where  t  and  d  are  the  coaponents  ofT  and  D,  respectively. 


Just  as  di 


i^erential 


foras  such  as  egaation  (33)  suggest 
integral  forms  like  equation  (35)  ,  corotational  derivatives  such 
as  equation  (46)  can  be  utilized  in  aaterially  objective  integral 
equations.  Such  equations  can  be  generalized,  in  a  Banner 
similar  to  the  generalization  of  the  corotational  Jeffreys  model 
to  the  Oldroyd  eight-constant  aodel,  to  give  the  Goddard 
memory-integral  expansion 


t  .  t  t  .... 

T  =  -  /  G^(t-t')C,dt*  -(1/2)/  /  G  (t-t», t-t")  [£*  .£"  ♦  £".£'  ]dt"dt* 

-  (1/2)/  /  /{2Giii  (t-t*  .t-t",  t-t  •)£•£••:£•  (49) 

♦  G^(t-t»,  t-t",  t-f)[£».  £".£•  ♦  £•.£".£•  ]}dt*dt"dt*  -  ... 

where  £"  is  a  generalized  strain  tensor  at  t",  and  the  G  ,  G  , 

G  ,  etc.  are  an  infinite  series  of  kernel  functions. 

Ill 

These  are  only  a  few  of  the  many  constitutive  equations 
which  have  been  foraulated.  Bore  extensive  compilations  can  be 
found  in  sources  such  as  Bird,  et  al^  [8].  Several  different 
aodels  are  still  necessary  because  no  single  formalization  has 
been  conclusively  deaonstrated  to  be  preeminent;  the  usefulness 
and  applicability  of  each  is  often  flow-specific  and  imperfect. 
Despite  these  limitations,  in  a  given  flow  situation  it  is  often 
possible  to  characterize  a  fluid  sufficiently  so  that 
semi-quantitative  agreement  and  prediction  can  be  made. 

Critical  evaluation  of  constitutive  equations  has  been 
hampered  considerably  by  experimental  constraints.  The  most 
common  experimentally  achievable  flows  are  vicometric.  These  are 
usually  approximations  to  steady  laminar  shear  flow,  with  the 
resultant  simplification  of  the  measurable  extra  stress  tensor  T, 
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see  equation  (39) .  The  degeneracy  in  £  results  in  a  degeneracy 
in  constitutive  models;  several  predict  the  same  behavior,  which 
leaves  no  basis  for  comparison.  This  has  led  experimentalists  to 
attempt  to  generate  non- viscometric  flows,  particularly 
elongational  flows.  A  uniaxial  elongational  flow  is  defined  by 
velocity  vector  and  rate  of  deforsation  tensor  of  the  fora 

I  =  [v^x  )/2,-v  (x  )/2]  (501) 

fi  -  ,  |  0  -1  0  |  <50B, 


where  «  is  the  elongation  rate.  One  characterization  of  a 
material  under  this  type  of  flow  is  the  elongational  viscosity  n 

t{11)  -  t  (22) 


n  - 
e 


e 
(51) 


The  value  of  the  elongational  viscosity  for  a  Newtonian  liquid  is 

three  times  the  value  of  the  shear  viscosity.  Rost  experimental 

techniques  to  measure  «|  have  been  restricted  by  mechanical 

e 

difficulties  to  the  testing  of  high  viscosity  polymer  melts  with 

6 

high  shear  viscosities,  »»  >10  poise,  at  low  rates  of  strain, 
.  -1  -1  o 

e<  10  sec  [22].  These  restrictions  do  not  apply  to  a  technique 
which  employs  bubble  dynamics  in  a  non-Ne wtonian  fluid. 

Pearson  and  Riddleman  [59,  60]  have  experimentally  observed, 
and  theoretically  analyzed,  the  collapse  of  a  spherically 
symmetric  cavity  in  a  viscoelastic  fluid.  In  spherical 

coordinates  this  is  easily  seen  to  be  a  uniaxial  elongational 

flow  and  kinematic  considerations  lead  to 


2R*R 


c  -  - 


(52) 


in  equation  (50b),  where  the  notation  of  section  II  is  retained. 
The  Rayleigh  equation  (4)  for  a  general  fluid  becomes 


P (8*8  ♦  ~-R2)  «  P(R>  -  P  ♦  /  [T.T]  dr  (53) 
2  -  R  r 

Various  mathematical  representations  of  £  can  then  be  used  to 
predict  bubble  behavior.  These  workers  chose  to  analyze  several 
models,  including  three-dimensional  versions  of  the  Haxwell 
equation  (32)  modified  with  codef ormational  (CD)  and  corotational 
(CR)  derivatives.  They  attempted  to  approximate  n  and  presented 
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(5«) 


resalts  in  terns  of  a  "bobble  pressure  function,"  e,  given  by 


•  *  "(3/2)  (p 


which  is  also  expressible  in 
teras  of  initial  values  of 
t(11),  t(22)  and  the  material 
constants  x  and  it  in  (32) . 
Experimentally,  the  apparatus 
schematically  represented  in 
figure  11  was  utilized  and 
bubble  behavior  recorded 
photographically.  The  results 
for  these  models  and 
experiaental  trials  are  shown 
graphically  in  figure  12  for  a 
solution  of  2.0}C  by  weight 
hydroxypropylcellulose  in 
water. 


I 


Stock  of  typical  tipcrimofital  bubbla  shop*  shoving  di 
lions  nsodad  for  tfco  data  analysis. _ 


FIGURE  11  [59,  P.718] 

As  e  increases,  the  agreement  between  theory  and  experiment 
progressively  worsens.  The  codef or national  nodel  fails  in  even  a 
gross  gualitative  sense  for  the  higher  elongational  rates.  This 
led  to  an  analysis  eaploying  wore  couplex  aodels,  including  a 
modified  corotational  Maxwell  fluid  and  an  integral  aodel,  which 
gave  much  more  accurate  predictions. 
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"3 


FIG.  12  1601 


♦(f)  data  compered  to  CD  (- - ) 

and  CR  ( - )  models  for  bubble 

collapse  in  a  HPC  solution,  t  m  0.28  »' 1 


♦(f)  data  compared  to  CD  ( - ) 

and  CR  ( - — )  models  for  bubble 

collapse  in  a  HPC  solution,  e  =  038  s“ 1 


♦(f)  data  compared  to  CD  ( - ) 

and  CR  ( - )  models  for  bubble 

collapse  in  a  HPC  solution,  i  m  0.77  s" 1 


♦(f)  data  compared  to  CD  (-  - — ) 

and  CR  ( — - )  models  for  bubble 

collapse  in  a  HPC  solution,  i  m  1.70s  1 


Won- Newtonian  effects  in  a  situation  in  which  departures 
froa  sphericity  becoae  important  have  been  analyzed  by  Vagner  and 
Slattery  [91]  as  an  extension  of  previous  work  for  Newtonian 
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droplets  in  a  uniform  Newtonian  flow  [79].  This  particular 
situation  involved  slow  flow  of  the  fluid,  which  allows  a 
perturbation  technique  to  include  inertial  effects.  A  second 
perturbation  was  then  utilized  to  include  non- Newtonian  effects. 
The  constitutive  sodel  employed  for  both  inner  and  outer  fluids 
was  the  Bivlin- Erick  sen  fluid  (43)  of  grade  3,  defined  by  (55) 


♦  1  ♦  a  A  *  4  k  *  4  k  *  4  tk 

1*1  2  2  3  1  4  3  5l*1 


♦  A  .4  ]  ♦  ♦  tr[A  ]4 
2  1  6  l*1  ia1 
(55) 


which  is  the  sinplest  grade  to  exhibit  both  noraal  stresses  and 
shear- dependent  viscosity.  Using  matched  asymptotic  expansions, 
drop  shape  was  found  to  be  calculable  through  numerous,  involved 
perturbation  coefficients.  These  shapes  progress  fron  "spherical 
to  prolate  spheroidal,  to  ovate  with  large  end  leading  and 
finally  to  a  teardrop  shape  with  extended  rear"  as  velocity 
increases,  figure  13.  These  shapes  are  in  qualitative  agreeaent 
with  experimental  observations.  Furthermore,  bubbles  have  been 
shown  to  be  a  special  case  included  in  this  analysis  [39]. 


Predicted  droptvt  ihopva. 


FIGURE  13  [91,  P.  1206] 


Even  an  accurate  constitutive  equation,  verified  by 
experiaent  for  aany  flows,  with  carefully  evaluated  constants,  is 
of  no  value  if  the  flow  geometry  of  interest  is  too  complicated 
to  allow  fruitful  analysis. 
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Such  is  the  case 
with  the  flows  which 
have  demonstrated 
the  non-lewtonian 
effect  of  immediate 
interest,  cavitation 
inhibition.  The 
results  of  Ellis,  gt 
ftl»  [ 27  ]  are  shown 
in  figure  1«  for 
flow  past  a 
hemispherical  nose, 
for  see  equation 
(1) .  Cavitation 
inception  was 
seasured  by 
detection  of  the 
initial  scattering 
of  a  laser  which  was 
adjusted  to  gra2ing 
incidence  on  the 
body.  The  effect 
was  noticeable  for 
concentrations  of 
polyethylene  oxide 
as  low  as  20  ppm,  as 
shown  in  figure  15. 


CONCENTRATION  OF  POLYETHYLENE  OXIDE  (W5R  SOM 
IN  RATER  IN  PARTS  PER  MILLION  SY  WEIGHT 


The  conditions  under  which  the  incipient  cavitation  number 
is  strongly  reduced  result  in  another  alteration  in  cavitation 
behavior  which  is  as  crucial  to  the  coaplete  understanding  of 
cavitation  inhibition.  The  appearance  of  the  subsequent  cavity 
is  changed;  the  void  appears  saaller  and  clearer  for  the  polyser 
solution  under  coaparable  flow  conditions. 
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This  change  is  noticeable  for  the 
flow  past  a  hemispherical  nose  used 
to  generate  the  data  of  figure  14, 
but  is  even  sore  apparent  for  a 
different  geometry.  In  figure  16 
the  voids  created  by  a  cylindrical 
cavitation  inducer  on  a  rotating 
disc  are  schematically  represented 
for  water  (a)  and  a  dilute  polyner 
solution  (b)  . 

Ting  [83]  takes  the  reduced  length 
of  the  inhibited  void  to  be  an  indication  of  a  "lower  cavitation 
intensity."  He  also  notes: 

In  water,  the  appearance  of  the  cavitation  bubbles  is 
very  violent  and  chaotic,  consisting  of  many  very  small 
bubbles.  As  a  result  it  scattered  much  of  the 
incipient  light  ...  However,  as  the  polyaer  is  added, 
the  cavity  looks  more  transparent,  and  shows  a  regular, 
saooth  wavy  pattern  at  the  vapor-liquid  interface  [83, 


V 


FIGORE  16  [83,  P.90C] 


P.900  ] 

This  indicates  that  a  polyaer  effect  is  not  liaited  to  the 
initial  stages  of  cavitation,  and  does  not  aerely  alter  a 
mechanism  of  nucleation.  A  strong  possibility  exists  that  a 
single  bubble  should  exhibit  coaaensurate  effects. 

The  simplest  model  system  which  might  exhibit  evidence  of 
cavitation  inhibition,  and  in  the  process  give  clues  to  the 
causes  of  the  phenomenon,  is  a  single  bubble  undergoing 
spherically  symmetric  dynamics  in  an  otherwise  quiescent  fluid. 
Pogler  and  Goddard  [29]  aay  have  been  the  first  to  analyze  this 
situation  for  a  non-Newtonian  fluid.  Their  analysis  began  with 
equation  (53),  which  can  be  transformed  to 


BP  ♦ 


3  .2 
— H 
2 


Pl-P- 


2  a 


- 


*  t(rr) 


dr 


p  pB  Hr 

A  siaple,  linear  viscoelastic  constitutive  equation  was  chosen 


(56) 


t  (t) 
rr 


t 

-2f  N  (t- 1  *)  d  ( t  •)  dt* 
o  rr 


(57) 
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and  the  aeaory  function  combined  viscous  and  "Haxwellian" 
coaponents  as 

*(t)  *  M«(t)  ♦  G  i-t/K)  (58) 

o 

Diaensionless  variables  and  a  characteristic  collapse  tiae  were 
defined  as 


t  =  a/R 


I  =  I(t  ) 
1  1 


t*  =  t/t 


t  =  R  (p/p  ) 
c  o 


1/2 


(59) 


with  the  equation  of  notion  becoaing 
P4.-P-  2  «Y 


YY  ♦ 


3  .2 
—  Y 

2 


P 

12HC,  t* 

—  t 


(exp[- 


1  v*«‘ 

(t»-t,  ) 


]J 


Y  Y  ln(Y1/Y) 


(60) 


*7- 


-dt 


1 


for  initial  conditions  Y(0)=1  and  Y(0)=0. 
parameters  arise  in  this  equation 

A 

I  =  - , 

De  t« 

G.t« 


Four  diaensionless 


a  Deborah  number 


N 


El 


Fe 


We 


an  elastic  nuaber 


a  Reynolds  nuaber 


a  Weber  nuaber 


(61) 


long  relaxation 
demonstrated  a 
large  Reynolds 


By  initially  focusing  attention  on  fluids  with 

times,  corresponding  to  V  the  authors 

Ce 

definite  elastic  effect  on  bubble  collapse.  For 

number  the  fluid  approaches  the  limiting  case  of  a  purely  elastic 

aaterial  and  rebound  short  of  collapse  is,  at  least 

theoretically,  possible.  Smaller  values  of  li  introduced 

R  e 

viscous  daaping  and  this  new  mode  of  energy  dissipation  increased 
the  value  of  the  radius  at  which  rebound  occurred. 

This  work  also  included  some  calculations  for  finite  Deborah 
nusbers.  In  this  case  the  void  will  eventually  collapse  to  zero 
radius.  However,  this  say  not  be  a  aonotonic  decrease  in  radios 
to  collapse  since  damped  oscillations  may  appear  superposed  on 
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the  overall  redaction  of  babble  size,  A  value  of  M  *0.51  is  the 

Oe 

lover  lisit  of  this  oscillatory  phenomenon  for  specified 

parameter  values  [I  *•,  p  /G  =10/7.  M  *•].  Thus  elasticity 

Re  •  o  He 

effects  can  have  a  pronounced  influence  on  the  behavior  of  a 
spherical  cavity  for  A/t*0(1). 

The  full  explanation  for  cavitation  inhibition  is  not  this 
simple  because  the  ratio  of  elastic  relaxation  time  to 
characteristic  collapse  tine  in  a  dilute  polymer  solution  is 
expected  to  be  such  less  than  1.  Ting  [82]  perforaed  a  aore 
realistic  analysis  by  employing  a  materially  objective 
constitutive  model,  the  Oldroyd  viscoelastic  fluid 


D  D 

t  ♦  A  - t  =  2ij[d.  .  +  a  - d  _  ] 

ij  1  Dt  ij  ij  2  Dt  ij 

0  3  3 

- b  =  —— b  ♦  v  — —  b  ♦  w  h  +wb  -  d  b  -d 

Dt  ij  3t  ij  k  3xfc  ij  it  kj  jk  ki  ik  kj  jk 


For  a  dilute  solution,  the  material  constants  in  equation 
can  be  related  to  aolecular  and  conponent  parameters  as 

n  =  n  (!  ♦  c[n]) 


(62A) 


b 

ki 

<62B) 

<62A) 


(63A) 


o 


A  *  A 
1 


(63B) 


(63C) 


A 

a2  ■  ~;?;y 

c[»j]  <  8  (6 3D) 

where  «i  is  the  solution  viscosity,  rj  is  the  solvent  viscosity,  a 

.  o 

is  the  terainal  relaxation  tine  of  polyaer  aolecules,  c  is  the 

polymer  concentration,  and  [»)]  is  the  polymer  intrinsic 

viscosity.  The  mechanism  by  vhich  bubble  growth  is  triggered 

fron  the  initial  eguilibriun  radius  R  is  a  step  change  reduction 

o 

in  the  ambient  pressure  of  magnitude  p*.  A  characteristic  time 

different  from  t  defined  by  Fogler  and  Goddard  (59)  is  employed 

c 

and  theraal  effects,  as  analyzed  by  Plesset  and  Zwick  [67],  are 
included.  The  new  dimensionless  variables  and  characteristic 
time  scale  are  given  by 
R 

~R  ~ 

*  (60) 


S  = 


an 


Pi(T.)-p^(0) 
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t 

t* 


t*  =  B  (p/p*) 
o 


1/2 


and  the  analogne  to  equation  (60)  becomes 
4 


••  3*2 

SS  ♦  — - S  ♦ 

2 


S  -3y  1  -1  1-3y 

-  =  an(1  -  S  ) - S  (1  -  S  )  ♦  H  (a) 

Re  S  8 

a  S^x)  S  (X)  a  S(X)  S*  (X) 

-  «  f  — — = — t - dx  -  2B  f  exp[  (x-at)  t'/A  ]  - [  1  ♦  — s - ]dx 

0  [J*S4(y)dy]  6  S(«)  S*(«) 

where  H  (a)  is  the  Heaviside  step  function,  y  is  the  polytropic 
gas  exponent,  and  the  new  dimensionless  parameters  are 


effects 


w  = 


*  = 


Re  = 


2  a 

aL p ' 


a  Reber  nuaber 

which  aeasures  theraal 


D  1/2  3/2 

( - )  t* 

w 


(P  *p) 


1/2 


(66) 


a  Reynolds  nuaber 


E  =  ctt> 


an  elastic  nuaber 


a  P* 

Numerical  calculations  of  R/R  vs.  tiae  were  carried  out  for 

parameter  values  corresponding  to  a  500  ppa  solution  of 

5 

polyethylene  oxide  (mw=4.5x10  )  for  initial  bubble  size 
1.02x10  cm  and  initial  temperature  T  =103  C. 


(65) 
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Besults,  figure  17,  show 

little  elastic  effect. 

Experimental  results  for  a 

similar  situation  [81] 

demonstrate  a  comparably 

small  effect.  The  gross 

cavitation  inhibition  present 

in  flow  systems  is  not  in 

evidence  here.  The  values  of 

H  used  by  Fogler  and 
De 

Goddard,  which  did  generate  a 
large  effect,  are  seen  to  be 
physically  unrealistic. 


ft/RCsse*1) 


riimt- 

2»ick 

Thsory 


V=o.* 


- -  Solvent 

-  Polyasr 


t(ssc)  - 


10-6  10-5  10-4 

.  Initial  bubble  growth  rate*  In  the  solvent  and  In  the 
polymer  solution. 

FIGURE  17  [ 82,  P.1430] 

Since  significant  cavitation  inhibition  is  not  present  in  a 
spherically  syaaetric  flow,  generation  and  analysis  of  aore 
complicated  bubble  geometries  is  indicated. 

The  jet-foraing  behavior  _ «1 _ 


of  a  collapsing  bubble  j  f  R  ' 

"* ;  1  - - - — •  Definition  of  sals 

near  a  solid  boundary  was  i — r  g  J  geometric  characteristics  of  s 

.  .  .  V  y  bubble  near  a  solid  wall. 

observed  by  Chahxne  and  i  !  N. 

*  i  ■ — 

Fruman  [14].  - 

FIGURE  18  [  14,  P.  1406  ] 
They  defined  geometric  characteristics  as  in  figure  18,  and  a 
ratio  t] 

ti  =  R  /L  (67) 

C,max 

The  addition  of  polyaer  solute  seemed  to  stabilize  the  sphericity 
of  the  cavity  and  to  reduce  the  intensity  of  the  re-entering  jet 
in  experimental  trials,  see  figure  14. 
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Thus,  although  the 
mechanism  is  far  from 
clear,  it  does  sees  that 
this  non-Newtonian  effect 
can  be  observed  in  single 
bubbles . 


Comparison  of  bubble  behavior  in  water  and  Polyox 
solution  for  equivalent  values  of  >). 


FIGURE  19  [ 19,  P. 1407] 
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If  PBELIHIHART  ANALYSES 


In  any  attempt  to  find  an  explanation  for  and  clarify  the 
aajor  factors  contributing  to  a  phenomenon  as  intricate  as 
cavitation  inhibition,  it  is  necessary  to  obtain  preliminary 
results.  These  results  are  expected  to  be  first  steps  and 
half-steps  toward  comprehension  of  the  techniques  necessary  for  a 
ccaplete  analysis.  They  often  ignore  elements  of  the  physical 
situation  which  are  integral  to  the  effect  of  interest  or 
virtually  duplicate  previous  work .  The  hope  remains  that  a  basic 
understanding  of  the  importance  of  a  given  factor,  by  its 
inclusion  or  exclusion,  will  ensue. 

The  complete  model  system  for  this  analysis  consists  of  a 
cavitation  bubble  changing  size  in  a  non-Newtonian  fluid,  while 
this  fluid  experiences  an  imposed  flow.  Preliainary  analyses  in 
spherical  bubble  dynamics  seek  confirmation  of  previous  works  on 
the  contribution  of  viscous,  inertial,  therual,  surface  tension 
and  elastic  effects  and  clarification  of  the  role  of  each.  An 
imposed  flow  creates  nonspherical  bubbles,  and  artificial  means 
may  be  employed  to  study  the  effect  of  nonsphericity  alone.  The 
imposed  flow  alone  can  also  be  studied,  and  night  give  some  basis 
for  choice  of  a  useful  constitutive  equation.  After  each  element 
of  the  system  is  examined  individually,  it  is  presumed  that  they 
can  be  combined  in  a  workable  analysis. 


PRELIHIBABT  BESOLTS 


Viscosity  and  Inertia  in  Spherical  Babble  Growth 


Meglecting  the  thermal  and  elastic  terms  in  Ting's  analysis 
[89],  his  resalt  (65)  becoaes 

CM  (B)  (C)  (D)  (E)  (P) 

3  .2  4  S  -3v  1  -1  1-3v 

ss  ♦  — s  ♦ - -  an  (1  -  s  ) - s  (1  -  s  )  ♦  h  {«) 

2  Re  S  B 

(65') 

where  S  (a)  is  a  dimensionless  babble  radius  which  is  a  function 
of  dimensionless  time  a.  The  first  two  terms  of  this  expression, 
(A)  and  (B)  ,  embody  inertial  effects,  term  (C)  includes  viscous 
effects,  tern  (D)  includes  the  internal  and  anbient  pressure 
difference,  tera  (E)  takes  surface  tension  into  account  and  the 
last  term,  (F)  arises  from  the  step  change  decrease  in  ambient 
pressure  of  nagnitude  p*,  which  is  assuaed  to  trigger  bubble 
growth.  (See  section  III  for  a  wore  complete  explanation.) 

Ignoring  the  inertial  terms,  (A)  and  (B)  ,  the  creeping  flow 
expression  for  radius  vs.  time  results.  In  real  tine,  for  all 
t  >  0, 

4 1'  1/B  13  1  -1 

t  - - - -  /  (x[  (-an  ♦  - )x  ♦  ( - )x  ♦  (1.0  ♦  an)])  dx 

Re  1/R_  W  W 

*  (68A) 

X  =  1/R  (68  B) 

There  exists  a  critical  value  of  the  initial  radius,  F  ,  which 

o 

resalts  when 


1 

- —  =  an  (69  A) 


C  2a 

F  =  - 

Physically,  this  is  the  eguilibriua  radius  for 
containing  no  noncondensable  gas;  i.e.  it  is 
physically  realizable  equilibrium  radius  for  a 
pressure,  anbient  pressure  and  surface  tension, 
equation  (68A)  can  be  simplified  to  qive  the  time 
growth  as 


a 

(69B) 

cavity 

the 

minimum 

given 

i  vapor 

For 

c 

R  =  R 

for 

0  o 

bubble 

t 


t» 

Be  (1*an) 


ln[ 


s  -  an/(i*an) 
1  -  an/(i*an) 


>3 


(70) 


c 

A  more  complicated  general  analytic  eolation  exists  for  R  >  R  . 

.  o  o 

Calculated  profiles  are  shown  in  figures  20 ,  21,  and  22  for  the 

o 

physical  parameters  corresponding  to  water  superheated  to  103  C. 

Retaining  tens  (A)  and  (B) ,  but  assuming  tera  (C)  is 
negligible,  the  result  for  the  inviscid  case  is 


t  = 


B  (fi/[p.  (T  )  -  p  10)  ]} 
o  l  o  *• 


1/2 


1/R  2  Re 4 

/  (—(1*1)  -  x  ♦  [ - 

1/R.  3  L  16W 


T*n 


-1/2  dx 


(71) 


c 

for  R  * R  .  This  is  an  elliptic  integral.  However,  analytical 
o  o 

results  reguire  numerous  tedious  transf oraations  to  arrive  at 
standard  forms,  and  tabular  results  of  the  subsequently  required 
values  for  the  three  types  of  elliptic  integrals  are  not  readily 
available.  Numerical  integration  was  employed  as  a  simple 
alternative  to  generate  the  profiles  of  figures  20,  21  and  22. 

Dergar abedian 's  experimental  results,  which  are  expected  to 
include  thermal  effects  and  which  closely  approximate  the 
analytic  results  of  Ting  (recall  Ting's  analysis  includes  an 
elastic  term) ,  see  figure  17,  are  also  shown  in  figures  20  and 
21.  From  the  lack  of  agreement  between  these  data  and  the 
theoretically  generated  results,  it  is  apparent  that  the  excluded 
effects  are  very  important.  Inertia  also  appears  to  be  more 
important  than  viscosity  in  this  situation. 
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EXPERIMENT 


Imposed  Plow 


To  investigate  the  effect  of  an  imposed  elongational  flow  on 
a  spherical  cavity,  an  attempt  was  made  to  calculate  the  initial 
deforaation  of  this  cavity  in  an  elastic  solid  and  also  in  a 
viscoelastic  liguid.  The  cases  considered  were  a  cavity  with 
constant  volume,  and  one  with  a  volume  profile  given  by  V  (t) . 
The  interior  fluid  was  assuaed  to  be  inviscid,  both  were 
incompressible,  and  inertial  effects  were  neglected  in  the 
eguations  of  aotion.  This  analysis  follows  the  work  of  Frohlich 
and  Sack  [32]. 

The  elongational  flow  is  assumed  to  be  undisturbed  at  r=R, 


where  R»1 


u  (R, 0,  t)  =  2e  (t)  RP  (cos  a) 
r  2 

u  (R, e,t)  =  e (t)  RP' (cos  a) 

0  2 


<72A) 

<72tf> 


in  spherical  coordinates,  where  P  is  the  second  Legendre 
polynomial 

The  foraalisa  of  Happel  and  Brenner  [38]  can  be  used  to 
generate  the  velocity  field  for  a  Newtonian  fluid  in  teras  ox 
spherical  haraonics.  The  aatheaatical  property  of  Legendre 
polynomials 


£  P .  («)  sin  a  da  *  0, 


aeans  that  no  net  change  in  cavity  voluae  can  result  froa  any 

term  in  a  series  expansion  for  v  except  the  spherically 

syaaetric  P  tern.  Thus,  the  radial  velocity  at  the  gas-liquid 
o 

interface  is  initially  given  by 

u  (r=a)  =  a (0)  (7tt) 

r 

i(0)  =  P(0)/(4,a  )  (75) 

for  any  given  volume  profile. 

Solving  for  the  velocity  field,  temporarily  ignoring  the 
zero  tangential  stress  condition  at  the  surface,  but  requiring 
that 

u  (r*a)  *  0  (76) 

a 
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(771) 


the  velocity  profile  exterior  to  the  spherical  body  is 
a1.  3aC  5a* 

a  =  — r®  ♦  [ - r -  r  ♦  2r]  *  (cos  a) 

r  r*-  rT  r  2 

c  . 

a  .  d 

u  =  (r - T-)e  [-—P  (cos  a)  ] 

a  r*  da  2 

Calculating  the  tangential  stress  at  the  surface 

d 

d  (r=a)  =  r  [—-P  (cos  a)  ]  f  0 
re  da  2 


(77B) 


(78) 


Additional  velocity  coiponents,  which  decay  as  r  increases, 
can  be  superposed  onto  this  solution  of  the  linear  equations  to 
give  a  coabination  which  satisfies  the  desired  zero  stress 

condition.  The  result  becoaes 

a  % 

a  .  a 

u  =  ——a  ♦  2*P  (cos  a)[-  ♦  rl  (79A) 

r  r*  2  r  *• 

.  r  d 

u  =  v  (—)  (cos  »)  (79 B) 

a  2  da  2 

this  axially  symmetric  velocity  field  can  be  recast  as  a  streaa 
function  given  by 


2.  3  3  r  .  2 

*  =  [a  a  ♦  (a  -r  )— — sin  a]  cos  t>  (80) 

2 

In  dimensionless  terms  this  becomes 

3  2 

*  =  *j[«  ♦  O-x  )  (1-m  )  ]  (80*) 

where 


v  -  cos  a 
x  =  r/a 
a  -  a/(ac) 
3. 


(81) 


a  =  */(a  r) 

Streamlines  for  various  values  of  a  are  shown  in  figures  23, 
24,  25,  6  26.  The  simplicity  of  this  solution  suggests  that,  at 
least  for  creeping  flows,  initial  probleas  are  tenable 
analytically.  Also,  the  form  of  the  solution  and  the  streaalines 
confirm  that  the  bubble  will  deform,  and  that  technigues  for 
nonspherical  surfaces  are  necessary  for  subsequent  work. 


A  Perturbed,  Growing  Cavity 


The  effect  of  nonsphericity  is  studied  through  the 
specification  of  an  artificial  situation.  Assume  an  inviscid 
cavity  is  generated,  with  a  fixed,  specified  shape  which  is  a 
slight  perturbation  from  a  sphere.  Also  assume  that  the  value  of 
the  eguivalent  radius  of  the  cavity  is  given  as  a  function  of 
time,  and  that  inertial  effects  are  negligible  in  the  Bewtonian 
surrounding  fluid. 

Specify  the  shape,  for  small  paraaeter  e  , as 

P  (t)  =  a(t)[  1  ♦  ,f  (M>  ]  (82 A) 


max[f  (M)  ]  =  0(1) 

u  =  COS  0 

tj  =  fluid  viscosity 
Perform  a  perturbation  analysis 


(0) 

(1) 
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(0) 
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p  =  p 

♦  r  p 

2 

♦  0(r  ) 


2 

♦  0(,  ) 


(82B) 

(82C) 


(83A) 

(83B) 


The  geosetric  perturbation  function  f  (»*)  can  be  expanded  as 

a  summation  of  surface  harmonics,  P  (m) •  If  a  trial  function 

n 

f  (u)  is  chosen 

=  P  2<*i)  (84) 

and  the  expansions  of  Cox  [20]  are  applied,  then  it  is  easy  to 
show  that 

a1  .  2  2  2 

a  =  — —a  ♦  r  (3aP  (  m)  (a/r)  [3(a/r)  -  2]}  ♦  0(r  )  (85A) 

r  rx  2 

•  a4  «/i  2 

m  *  0  ♦  r[  3a — — PMu-n  )  ]  ♦  0(r  )  (85B) 

0  T*  2 

The  homogeneous  nature  of  the  boundary  condition  applied  at  the 
surface,  zero  tangential  stress,  means  that  this  velocity  profile 
is  also  valid  for  a  linear  viscoelastic  fluid.  Again,  as  in  the 
previous  result,  the  creeping  flow  problea  has  been  shown  to  be 
straightforward  and  amenable  to  a  perturbation  analysis. 
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PROPOSED  PBELIHI BABY  1HALYS ES 


Elasticity,  Surface  Tension  and  Theraal  Effects 

The  result  of  Pogler  and  Goddard,  equation  (60) ,  and  of 
Ting,  equation  (65) ,  can  be  solved  for  more  parameter  values  to 
generate  profiles  like  those  in  figures  20,  21  and  22.  With 
carefully  chosen  values  for  these  physical  constants,  i.e. 
systematic  departures  from  the  values  expected  for  pure  water, 
dilate  polyaer  effects  on  spherically  syaaetric  growth  are 
expected  to  be  denonstrated  to  be  ainiaal. 

Imposed  Flow 

First,  the  analyses  in  the  prelininary  results  section  under 
this  subheading  require  coapletion.  In  addition,  the  initial 
deformation  of  a  spherical  cavity  under  other  imposed  flows,  such 
as  simple  shear  or  the  combination  of  shear  and  vorticity 
afforded  by  an  orthogonal  rheoneter,  can  be  calculated.  The 
stress  field  should  also  be  evaluated  with  particular  attention 
paid  to  changes  in  the  magnitude  of  noraal  stresses.  These 
attenpts  should  eaploy  various  constitutive  relations,  with 
parameters  corresponding  to  liquids  ranging  from  a  purely 
viscous,  to  a  coapletely  elastic  aaterial. 


f  ULTIMATE  PLUS  and  GOALS  for  this  Besearch 

The  preferred  goal  of  this  research  is  to  determine,  by 
aeans  of  theory  and  experiment,  whether  a  aodel  systea  consisting 
of  a  single  babble  in  a  well-characterized  non-Newtonian  flow 
exhibits  the  gross  effects  of  cavitation  inhibition. 

THEORETICAL 

Host  of  the  preceding  pages  have  been  devoted  to  an 
introduction  to  the  approaches  and  -.echnigues  available  for 
theoretical  analysis  of  a  bubble  system.  At  this  stage,  it  is 
iapossible  to  predict  which  single  technigue  will  be  most 
important.  All  will  probably  become  necessary. 

However  it  is  achieved,  the  resultant  theoretical  system 
should  be  flexible  and  capable  of  including  various  flow 
configurations.  This  is  expected  to  allow  some  evaluation  of  the 
importance  of  various  factors  within  a  flow,  such  as  vorticity 
and  shear  rate.  A  procedure  which  nay  be  particularly  important 
is  a  search  for  an  imposed  flow  which  allows  a  steady  state 
nonspherical  shape  of  a  bubble  with  constant  internal  pressure. 
The  existence  of  such  a  flow  is  suggested  by  the  first  special 
case  mentioned  by  Ballison  [71],  the  "weak  flow"  conditions.  For 
a  bubble,  the  viscosity  ratio  a  (23a)  is  very  small,  so  the 
reguisite  flow  strength  for  eguilibrium  might  be  too  low  to 
generate  important  nonspherical  effects.  If  this  is  not  the 
case,  it  is  possible  that  the  change  in  cavitation  intensity 
noted  by  Ting  [82]  will  be  manifest  in  a  change  of  share  and/or 
volume  of  the  steady  state  cavity.  This  is  likely  to  be  a  mucn 
simpler  analysis  than  any  including  bubble  growth,  especially  if 
ellipsoidal  harmonics  are  applicable. 

Just  as  Pearson  and  Middleman's  analysis  gives  some  basis 
for  the  evaluatiob  of  constitutive  relations,  the  results  here 
may  be  be  of  value  in  the  assessuent  of  constitutive  accuracy. 
Even  without  experimental  results  for  comparison,  physically 
unreasonable  results  may  eliminate  certain  eguations.  with 
experimental  data,  still  closer  scrutiny  should  be  possible. 

The  problems  expected  here  are  numerous,  particularly  since 
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the  nonlinearity  of  the  complete  equations  of  motion  makes 
analytic  solutions  unlikely.  Exact  solutions  for  limiting  cases 
will  be  found  whenever  possible,  and  the  need  for  a  great  deal  or 
effort  and  cleverness  is  anticipated  in  order  to  qenerate 
numerical  solutions. 

EXPERIMENTAL 

Two  major  types  of  experiments  are  envisioned  in  plans  for 
this  research:  exploratory  and  confirmatory.  Exploratory 
experimentation  is  conducted  concurrent  with  theoretical  work. 
Its  purpose  is  to  guide  the  development  of  theory,  eliminating 
unimportant  approaches  and  distinguishing  physical  effects  fron 
purely  mathematical  ones.  From  this  initial  type  of  experimental 
work  and  theoretical  analysis,  a  full  theoretical  model  should 
arise,  and  final  testing  and  adjustment  is  made  through 
"confirmatory"  experimental  trials. 

Before  any  of  this  work  begins,  logistical  and  technical 
problems  must  be  solved  or  circumvented.  These  include  choice  of 
materials  and  apparatus,  procurement  of  tunds,  and  even  ordering 
and  deliverv  constraints.  Once  the  apparatus  is  assemble  1, 
technigues  must  be  learned  and  noned,  and  trial  runs  conducted. 
A  long  process  is  expected  to  precede  any  experiments  whicn 
approach  the  model  system.  Initial  decisions  concerning 
apparatus  will  depend  on  preliminary  analyses  and  past  work  for 
guidelines.  Careful  choice  and  characterization  of  lluids  must 
te  made,  then  means  must  be  devised  for  flow  generation,  buDble 
generation,  and  the  recording  and  analysis  of  bubble  behavior. 

The  reference  fluid  will,  of  course,  be  water,  while  the 
test  fluids  are  expected  to  be  dilute  solutions  of  polymers  such 
as  polyethylene  oxide  [14,  27  J,  guar  gum  (.27],  or  polyacrylamide 
[  15,  83].  The  choice  will  depend  upon  the  behavior  of  candidate 
additives  under  some  new  conditions  which  may  be  present  in  the 
experimental  procedure.  All  fluids  must  be  characterized 
viscometrically ,  and  also  checked  for  variations  in  physical 
parameters,  such  as  surface  tension,  from  the  values  for  pure 
water.  The  flow  apparatus  used  to  qenerate  the  imposed  flow  may 
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also  be  useful  for  fluid  characterization 


The  specific  aachinery  necessary  for  flow  generation  will  of 
course  depend  on  the  desired  flow.  A  simple  one -dimensional 
shear  flow  can  be  generated  by  a  two-belt  apparatus  (see  figure 
27a).  A  two-dimensional  extensional  flow  results  front  a 
four-belt  apparatus  (see  figure  27b) .  A  Couette  flow  can  easily 
be  generated  experimentally  (see  figure  27c).  Perhaps  the  lost 
flexible,  single  flow  qeometry  is  present  in  an  orthoaonal 
rheometer  (figure  27d)  .  By  varying  the  offset,  rotation  rate  and 
plate  gap,  the  amount  of  vorticity  and  shear  present  in  the 
undisturbed  flow  can  be  varied  independently.  Hakimi  ana 
Schowalter  [37]  used  just  such  a  device  in  their  experiments  with 
drop  deformation.  Thus  tne  technology  for  and  components  of  a 
unit  should  be  availaole  in  the  Princeton  University  Department 
of  Chemical  Engineering,  where  this  research  is  being  conduct-i. 

(a)  (b)  (c)  (d) 

FIGURE  27  Possible  Experimental  Flow  Apparatus 

The  most  versatile  and  sophisticated  of  the  devices  which 
may  be  applicable  to  flow  generation,  and  will  definitely  be  the 
major  means  of  rheological  characterization,  is  a  FHEOflETRICS 
Mechanical  Spectrometer  ,  whicn  is  now  under  order  by  this 
department.  Several  flow  geometries  are  available  in  tnis 
device,  including  Couette  flow  and  orthogonal  rheometer  flow, 
with  a  wide  range  or  operating  parameters.  This  instrument  has 
the  added  advantage  of  capability  for  making  sensitive  dynamic 
measurements  of  macroscopic  elastic  effects  such  as  the  total 
normal  force  on  a  rheometer  plate.  By  comparing  these 
measurements  for  identical  imposed  flows,  with  ana  without  a 
cavity  present,  additional  conclusions  with  respect  to  concepts 
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such  as  flow  alteration  may  be  possible. 

Bubble  generation  can  take  on  two  levels  of  sophistication. 
If  a  steady  state  bubble  shape  is  found  to  be  achievable,  then 
cavity  formation  for  such  a  flow  may  be  as  simple  as  injection  of 
a  measured  volume  of  vapor  or  gas  into  the  static  fluid. 

To  produce  cavitation  bubbles  for  dynamic  study  the 
technique  with  the  largest  probability  of  success  is  the  use  of  a 
laser  to  trigger  nucleation.  A  Q-svitched  ruby  laser  has  been 
used  to  successfully  generate  volumes  in  pure  water  which  are 
well-modelled  as  single  cavitation  bubbles  [52,  53].  The 

0-switch  is  a  shutter  which  contracts  the  pulses  of  the  ruby 
laser  into  smaller  segments  of  length  30-50  nsec.  In  this  short 
time,  about  1  J  of  energy  is  introduced  into  a  small  region  of 
the  fluid.  Available  references  on  lasers  [51]  show  a  wide 
variety  of  similar  devices  to  be  evaluated. 

The  recording  system  is  also  very  different  for  the  steady 
state  or  dynamic  bubble.  The  steady  state  shape  requires  only 
sinqle  photoqraphs  of  two  or  three  views  of  the  bubble.  Hith 
adequate  lighting,  a  simple  Polaroid  camera  should  suffice, 
although  slide  capability  would  be  desired  to  allow  projection  of 
the  images. 

The  dynamic  bubbles  are  expected  to  change  on  a  time  scale 

of  about  10  sec.  or  less  (figures  20  and  21).  Recording  of 

such  rapid  events  photographically  requires  a  very  high  framing 

rate.  Cameras  with  rates  this  high,  and  higher,  are  available 

[75],  and  have  been  used  to  record  bubble  dynamics  (figure  28). 

These  cameras,  however,  are  extremely  expensive,  e.g.  the  Cordin 

6 

310  with  framing  rate  of  2x10  pictures  per  second  costs 
approximately  $100,000.  Rental  of  such  units  is  also  available. 
This  is  an  attractive  proposition,  especially  for  the  early 
trials  when  needs  are  not  completely  known.  If  the  flow  is  not 
axisymmetric,  more  than  one  view  of  a  bubble  will  be  necessary 
for  analysis.  Should  reproducibility  be  found  to  be  excellent, 
then  photography  of  the  different  views  in  different  trials  will 
te  acceptable;  if  not,  an  optical  system  for  recording  multiple 
images  on  one  frame  may  be  necessary. 
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Dvnsmics  of  a  lascr-produoed  spherical  bubble  near  a  eolid  boundary.  The 
framing  rate  "is  76000  frames/s,  the  maximum  bubble  radiua  R _ =  2-0  mm,  the  dis¬ 

tance  of  the  bubble  centre  from  the  boundary  £  =  4-0  mm  and  the  size  of  the  individual 
frames  is  7-2  x  4-6  mm. 


FIGURE  28  [ 53,  P.400] 

Analysis  of  the  massive  photographic  record  anticipated  will 
probably  require  countless  projections  of  individual  images,  and 
subsequent  fitting  to  analytical  expansions.  Fortunately, 
exploratory  trials  will  probably  be  most  useful  for  qualitative 
results,  reducing  the  need  for  exhaustive  examination  of  the 
images.  Later,  a  good  theoretical  model  will  give  accurate 
predictions,  minimizing  the  deviation  from  predicted  shapes  and 
the  need  for  fitting. 
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The  proposed  overall 
system  is  diagramaticall y 
presented  in  figure  29. 
Of  course,  the  details 
are  dependent  upon  future 
developments. 
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Block  diagram  of  experimental  setup. 

FIGURE  29  C  5^,  P.  2d ] 

Many  unforeseeable  problems  are  expected  to  arise  in  the 
course  of  this  work.  Others  may  be  at  least  partially 
anticipated  and  must  be  resolved  for  success  to  ce  achieved. 
Possible  problems  include: 

1)  Multiple  bubbles 

The  ruby  laser  is  thought  to  cause  nucleation  on  dust 
particles  [77].  The  polymer  solute  molecules  may  act  as 
nucleatincj  centers,  the  result  being  formation  of  multiple 
biinbles.  A  possible  solution  would  be  to  split  the  laser  ani  to 
focus  multiple  beams  on  a  single  point  from  many  angles. 

2)  Polymer  degradation 

The  laser  may  also  cause  thermal  degradation  cr 
polymer.  This  is  very  likely,  at  least  in  the  immediate  vicinity 
of  th®  focai  point,  and  must  be  evaluated  for  macroscopic  effect. 

3)  Surface  tension 

Impurities  in  water  are  thought  to  change  the  dynamics 
of  a  single  bubble  in  Stokes  flow  by  altering  surface  tension. 
Polymer  molecules  might  cause  a  similar  phenomenon. 

4)  Vapor  pressure 

Is  the  vapor  pressure  of  a  dilute  solution  the  same  as 
that  of  the  solvent?  Measurements  need  to  be  taken  here.  one 
possibility  for  testinq  the  coupled  effects  of  3)  and  4)  is  a 
static  bubble  comparison  between  water  and  the  solution,  since 
surface  tension  and  vapor  pressure  are  the  major  determinants  of 


equilibria!  radius. 

5)  Thermal  effects 

Error  nay  be  introduced  into  the  system  by  a  change  in 
theraal  properties  of  the  liquid  or  laser  heating  effects.  Heat 
effects  from  photographic  lighting  can  be  minimized  using  proper 
filters. 

6)  Hall  effects 

By  varying  apparatus  or  bubble  dimensions,  it  should  be 
possible  to  estimate  the  effects  of  the  finite  extent  of  the 
experimental  fluid.  In  using  the  Mechanical  Spectrometer,  normal 
force  readings  nay  also  give  an  indication  of  the  dirterence 
between  theoretical  flows  imposed  at  infinity  and  real  flows. 
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